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Electric field noise is a hinderance to the assembly of large scale quantum computers based on 
entangled trapped ions. Apart from ubiquitous technical noise sources, experimental studies of 
trapped ion heating have revealed additional limiting contributions to this noise, originating from 
atomic processes on the electrode surfaces. In a recent work [A. Safavi-Naini et al, Phys. Rev. A 
84, 023412 (2011)] we described a microscopic model for this excess electric field noise, which points 
a way towards a more systematic understanding of surface adsorbates as progenitors of electric 
field jitter noise. Here, we address the impact of surface monolayer contamination on adsorbate 
induced noise processes. By using exact numerical calculations for H and N atomic monolayers on 
an Au(lll) surface representing opposite extremes of physisorption and chemisorption, we show 
that an additional monolayer can significantly affect the noise power spectrum and either enhance 
or suppress the resulting heating rates. 



I. INTRODUCTION 

Ion trap miniaturization and a precise control of errors 
in the entangled qubits are two key prerequisites for us- 
ing trapped laser-cooled atomic ions as multi-qubit logic 
gates in a scalable quantum architecture ^-[3] • One main 
source of qubit error in such systems is the motional jit- 
ter of the collective behavior of the ions in micro traps 
[4T-H5] . Some of the unwanted heating noise is naturally 
mitigated by operating the traps at cryogenic tempera- 
tures [9j [TD] . However the noise still remains larger 
than the expected Johnson noise for the traps. Early 
observations of the dependance of the heating rate on 
position of the ions above the trap (~ 1/d 4 ) (I] and on 
the elapsed time in the ion- loading region [JJ |S] , lend cre- 
dence to the role played by surface contaminants. Further 
experiments with superconducting traps [10] corroborate 
the understanding that the noise source lies on the sur- 
face and not in the bulk. The confirmation has come 
more directly from two recent and complimentary exper- 
iments, where laser cleaning [12 and ion beam bombard- 
ment [13] of the trap electrodes led to a reduction of the 
noise. The experiment reported in |13j also identified the 
surface contaminants as carbon based, in the form of 2-3 
monolayers (MLs) of hydrocarbons. 

Theoretical studies of the anomalous heating have been 
largely phenomenological, aiming to explain the signa- 
tures of this noise. These models use the concept of patch 
potentials developed by Turchctte et al. to explain the 
motional heating in ions [H QTJ [T3J [T5] . In a recent work 
|16j we developed a microscopic model to predict the fea- 
tures of electric field noise (distance, frequency and tem- 
perature dependencies) from the details of atomic surface 
processes. This model is predicated on the idea that the 
noise in ion traps emanates from a random distribution 
of fluctuating dipoles associated with individual adatoms 
on a metallic electrode surface. 



In this work we extend our earlier theoretical treatment 
by investigating the dependence of this surface noise on 
the presence of an additional ML of atomic species on 
a gold surface. To do so we present detailed numeri- 
cal calculations on the adsorbate surface potential in the 
presence or absence of MLs with different reactivity. To 
this end, we chose He as the ML atom to represent weak 
binding (physisorbed species) and N as the ML atom 
for strong binding (chemisorbed species). For compu- 
tational simplicity we choose hydrogen (H) as our adsor- 
bate and compare the binding potential of H/Au(lll), 
H/He/Au(lll) and H/N/Au(lll). Using density func- 
tional theory (DFT) we obtain accurate data for these 
surface potentials, induced dipole moments and modifi- 
cations on the phonon density of states. Combined with 
the semi-analytical treatment of Ref . [TB] we use this data 
to extrapolate the resulting impact on phonon induced 
dipole fluctuations and find that MLs of different reactiv- 
ity can lead to the completely opposite effects of reducing 
or enhancing the noise. 

The remainder of the paper is structured as follows. 
In Sec. [H] we briefly review the problem of anomalous 
heating in ion traps and summarize the basic assumptions 
of the model detailed in Ref. [IB]. As the main part of 



this work we present in Sec. Ill our numerical results on 



adatom surface potentials and induced dipole moments 
for two different types of MLs on Au. Finally, in Sec. |IV| 
we discuss the impact of these findings on the adatom 
dipole fluctuation spectrum and conclude in Sec. [V] 



II. ANOMALOUS ION HEATING FROM 
FLUCTUATING SURFACE DIPOLES 

In micro-fabricated surface ion traps, which are cur- 
rently developed for quantum information processing, 
single or multiple ions are trapped by electric potentials 
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at a distance d of a few 100 /im above a metal electrode. 
The resulting trapping frequencies are typically around 
uj-t ~ 1 MHz and allow efficient laser cooling and coherent 
manipulations of the trapped ion. However, when cooled 
close to the quantum ground state of the trap, the ion 
motion is still associated with a comparably large elec- 
tric dipole moment dj w qag, where q is the charge and 
clq = y / 'h/(2mjujt) the zero-point motion for an ion mass 
to/. Fluctuating electric fields from the environment cou- 
ple to this dipole moment and excite the ion motion with 
a characteristic heating rate [I] 
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where Se(u) is the fluctuation spectrum of the electric 
field at the position of the ion. Since I\ limits the time 
for performing coherent manipulations of the ion, a de- 
tailed understanding of Se{w), its distance, frequency 
and temperature dependence, is of central importance for 
a further optimization and miniaturization of ion micro- 
traps. 

As the trapping distance d is decreased the ion be- 
comes increasingly more sensitive to electric noise emerg- 
ing from microscopic processes on the surface. In [16j we 
developed a microscopic model to describe the electric 
field noise, which is generated from a random distribu- 
tion of adatoms on a gold surface. In this case the field 
fluctuation spectrum for a planar trap geometry is given 
by 
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where a is the surface density of dipoles and Sn(w) — 
dr(/z z (r)/Lt z (0))e l " T is the spectrum of an individ- 
ual fluctuating dipole. Eq. |2]) predicts the expected d~ 4 
scaling and together with Eq. |l]) it relates the ion heat- 
ing rate to the microscopic dynamics of individual surface 
impurities. 



A. Phonon induced dipole fluctuations of adatoms 
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FIG. 1. (a) A schematic binding potential U[z) of an adsor- 
bate on a bare Au surface as a function of the adatom-surface 
distance z. The dotted lines are the bound states of the po- 
tential and the corresponding wavefunction is shown using 
solid lines, (b) The typical dependence of the phonon-induced 
dipole fluctuation spectrum of the adatom as a function of 
oj/Vq, where To is the characteristic transition rate from the 
first excited state to the ground state. The temperature is 
given in units of ^io, the separation between the ground and 
first excited vibrational states. See text and Ref. [16] for more 
details. 



Fig. [T] (a) shows a typical adatom-surface potential 
U(z), which is attractive at large distances z and has 
a sharp repulsive wall when the electronic wavefunctions 
of the adatom and the surface atoms start to overlap. 
The adatom-surface interaction is associated with a dis- 
tortion of the electronic wavefunctions which results in 
an induced dipole moment /-t z (z) perpendicular to the 
surface. At large distances one expects fi(z) ~ 1/z 4 [T7] 
and fj,(z ~ zq) can reach several Debye when the adatom 
touches the surface. 

The potential U (z) usually supports several bound vi- 
brational states | n) with vibrational frequencies v n and 
the adatom can undergo phonon induced transitions be- 
tween those vibrational states. For n > to the corre- 
sponding transitions rates are approximately given by 



r„ - ng{ " nm) \(n\U'(z)\m)\ 2 (n(u nm ) + 1) , (3) 
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Here M is the surface atom mass, g{oS) is the phonon 
density of states (PDOS) and n{u) = i/^/^bT) _ ^ 
the thermal phonon occupation number, which are both 
evaluated at the vibrational transition frequency v nm — 
v n — v m > 0. Due to the different average dipole moment 
fJ>n = ( n \^( z )\ n ) associated with each vibrational state, 
absorption and emission of phonons creates a fluctuat- 
ing dipole moment /z(f) as the adatom jumps between 
different levels In). 
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B. Dipole fluctuation spectrum 

From the above considerations and a detailed knowl- 
edge of the adatom-surface potential U(z), the induced 
dipole moment (i(z) and the phonon density of states 
g(u>), we can evaluate the dipole fluctuation spectrum 
and thereby the corresponding ion heating rate I\. The 
fluctuation spectrum is given by 

/CO 
dr((^(T) Mz (0))-(/i z (0)) 2 )e— , (5) 
-co 

where p z = ^ n p z ,nPn and p n — \n)(n\ is the projection 
operator on the vibrational level \n). 

The typical dependence of S'p(w) on frequency and 
temperature is shown in Fig. [ljb) establishing the rate 
r = Ti^ Q (T = 0) and the frequency v w = (Ei — E )/H 
as the relevant scales in the problem. In |16j we used ap- 
proximate analytic model to estimate the relevant scales 
for a broad range of adatom species, but assuming a clean 
gold surface. By using a harmonic approximation for 
U(z) and assuming that g(oj) ~ uj 2 we obtain [TB] 

^10~CW 2' r 0~^ 3 ' ( 6 ) 

where £ ~ O(l) is a numerical constant, m is the mass 
of the adsorbate, p is the density of the slab, Uo is the 
potential depth and v is the speed of sound. 

In the following we consider a more realistic scenario 
and evaluate the potential modifications of the dipole 
fluctuation spectrum due to the presence of an additional 
ML of atoms on top of the Au surface. To do so we 
present in the following section exact numerical calcula- 
tions adatom-surface potentials for the case of He and N 
monolayers, which provides us with an estimate for the 
minimal and maximal expected modification of the sur- 
face potential. Approximate analytic expressions for r 
given in Eq. |6]) allow us to extend these predictions for 
various monolayer-adatom combinations. 

III. ATOM-MONOLAYER-GOLD SURFACE 
INTERACTION 

The asymptotic potential for a polarizable atom with 
dynamic polarizability which approaches a surface 

of dielectric constant e is U(z ^ zq) ~ — ^j-, where 

C*3 — j^: / a(iuj)dw and z is the normal to the surface. As 
the atom approaches the surface, the interaction energy 
increases due to exchange of energy with phonons, lead- 
ing to adsorption at some equilibrium distance zq near 
the surface, beyond which, for shorter distances of ap- 
proach, the energy cost of electronic exchange between 
the electrons of the adsorbate atom and the bulk atom 
becomes too great to overcome, leading to a repulsive 
wall in the interaction potential. Below, we present ab 
initio calculations of the interaction potential energies 



normal to the substrate surface of H atoms with a ML of 
adsorbate atoms He and N on top of the Au(lll) surface. 



A. ML surface interaction potentials 

Total-energy calculations of bulk Au and Au(lll) sur- 
faces, with and without He and N adsorbate atoms, 
were performed using the spin-polarized density func- 
tional theory as implemented in the Vienna Ab initio 
Software Package (VASP) [18]. The exchange correlation 
energy was calculated using the local gradient approxi- 
mation (LDA) with the parametrization of Perdew and 
Wang (PWC) US]. 

The interaction between valence electrons and ionic 
cores was described by the Projector Augmented Wave 
(PAW) method [2JB HQ- The Au 5d 10 6s\ N 2s 2 2p 3 and 
He Is 2 electrons were treated explicitly as valence elec- 
trons in the Kohn-Sham (KS) equations and the remain- 
ing cores were represented by PAW pseudopotentials. 
The KS equations were solved using the blocked Davidson 
iterative matrix diagonalization scheme followed by the 
residual vector minimization method. The plane-wave 
cutoff energy for the electronic wavefunctions was set to 
500 eV. 

All structures were optimized with periodic bound- 
ary conditions applied using the conjugate gradient 
method, accelerated using the Methfessel-Paxton Fermi- 
level smearing [22] with a width of 0.2 eV. The total en- 
ergy of the system and Hcllmann-Feynman forces acting 
on atoms were calculated with convergence tolerances set 
to 10 -3 eV and 0.01 eV/A, respectively. Structural opti- 
mizations and properties calculations were carried out us- 
ing the Monkhorst-Pack special fc-point scheme [23] with 
llxllxll and 7x7x1 meshes for integrations in 
the Brillouin zone (BZ) of bulk and slab systems, respec- 
tively. 

The supercell consisted of a three-layer thick gold slab 
with (111) orientation and a p(2 x 2) mesh unit, cov- 
ered by He or N adsorbate atoms on one side of the 
slab model. The calculated lattice constant of bulk Au 
was 4.06 A, in close agreement with the experimental 
value of 4.0780 A at 25°C [M]. The lattice parameters 
of the p(2 x 2) Au(lll) surface constructed by cleaving 
the optimized bulk structure were a = b = 5.74 A and 
c = 25.00 A, with ca. 20.00 A vacuum separating slabs, 
and a = /? = 90° and 7 = 120°. Although a large vacuum 
region (ca. 20 A) was used between periodic slabs, the 
creation of dipoles upon adsorption of atoms on only one 
side of the slab can lead to spurious interactions between 
the dipoles of successive slabs. In order to circumvent 
this problem, a dipole correction was applied by means 
of a dipole layer placed in the vacuum region following 
the method outlined by Neugebauer and Schcfflcr [25] . 

The He/Au(lll) and N/Au(lll) interaction potentials 
were calculated by gradually moving a single He or N 
atom along the z-axis normal to the Au(lll) surface. 
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B. Atomic adsorption on Au(lll) 

Four different atomic adsorption sites are possible onto 
a Au(lll) surface: 1) a bridge site between two gold 
atoms, 2) on top of a gold atom, 3) in a hollow site 
between three gold atoms, termed an hexagonal close 
packed (hep) site when there is a gold atom in the layer 
directly beneath the surface layer, or 4) termed a face- 
centered cubic (fee) site when there is a hole in the layer 
directly beneath the surface layer. 

Total-energy calculations indicate that a single He 
atom adsorbs preferentially at the bridge site (E = 
—48.760 eV), slightly more energetically favorable than 
at the top site (-48.756 eV), the fee site (-48.755 eV) 
and the hep site (—48.749 eV). The elongated equilib- 
rium He-Au bond distance of 3.58 A suggests that He 
at the bridge site is weakly physisorbed to the Au(lll) 
surface. For the adsorption of a single N atom, the fee 
site is energetically preferred (E = —70.393 eV) over the 
hep site (E = -69.792 eV), the bridge site (-69.098 eV) 
and the top site (E — —66.510 eV). Contrasting with the 
He adsorbate, the N atom occupying the fee site appears 
chemisorbed to the Au(lll) surface with a short N-Au 
bond distance of 2.05 A. 

Fig.[2]shows the adsorbate potentials for the two mono- 
layers as well as the bare Au surface. Note that we have 
shifted the potentials so that E(zq — > oo) = 0. The 
presence of the weakly adsorbed He results in a much 
shallower potential while the chemisorbed N has the op- 
posite effect, creating a deeper well that support more 
bound states. 
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FIG. 2. The binding potentials for H adsorbate atoms on 
bare Au surface (dotted) and Au surface covered with He 
(solid) or N (dashed) monolayers. The peak of the local po- 
tential at the position of the first Au layer, in direct contact 
with the N(He) ML, does (not) vary appreciably indicating 
the formation of a sizable (negligible) dipole moment. 



C. Atomic diffusion on Au(lll) 

Although our primary focus in this work is on dipole 
fluctuations induced by atomic motion perpendicular to 
the surface, the numerical data obtained in the previous 
part allows us with no additional effort to evaluate the 
diffusion rates of adatoms parallel to the surface. Such 
data could be relevant to other, diffusion related noise 
processes [2"6"r[2"5] under identical conditions. 

The surface diffusion coefficient is given by D = ^-r^T 
for an fcc(lll) surface, where tq is the lattice parameter 
[2.87 A for Au(lll)] and T is the jump rate. Two differ- 
ent diffusion regimes exist, namely a thermally activated 
regime and a quantum tunneling regime. 

In the thermally activated regime, the thermal jump 
rate, r t h om , can be calculated through the Arrhenius 
formula [2"9"| . 

Tthorm = 7exp (-j^TpJ > ( 7 ) 

where E a is the activation energy barrier, fee is the Boltz- 
mann constant, T the temperature, and 7 is a prefactor 
which contains dynamical quantities; 7 rs 10 12 — 10 13 Hz 
for most surfaces [3U]- The diffusion of a He atom on 
Au(lll) is expected to be nearly barrierless due to their 
weak physisorption interaction. 

For the diffusion of a single N atom between adja- 
cent stable fee sites on Au(lll) the activation energy 
was calculated in this study to be E a = 0.17 eV using 
the nudge elastic band method within the DFT frame- 
work. This result is close to the value of ~ 0.10 eV 
determined previously from experiment and theory for 
the diffusion of Cr on Au(lll) surface [3T]. Following a 
simple Redhead's analysis of the migration temperature 
on solid surfaces [32], we estimate the temperature nec- 
essary for a N adatom to overcome this energy barrier to 
be « 65 K [E = 0.06 T kcal mol" 1 K" 1 ]. Let us note that 
Redhead's law predicts the migration temperature of Cr 
adatoms on Au(lll) to be 39 K, in excellent agreement 
with experimental findings |31j . The thermally-activated 
diffusion jump rate of N on Au(lll) at 65 K is calcu- 
lated to be r t hcrm = 0.3 Hz for an activation energy of 
0.17 eV and a prefactor 7 = 5 x 10 12 Hz; this corresponds 
to a surface diffusion coefficient D = 1.2 x 10~ 16 cm 2 s _1 . 
In the temperature range 60 — 70 K, r t hcrm and D vary 
from 2.6 x 10~ 2 to 2.9 Hz and from 9.3 x 10" 18 to 
1.0 x 10 -15 cm 2 s _1 , respectively. 

Below w 65 K, classical thermally activated surface 
diffusion can be excluded as the origin of the N adatom 
diffusion and the site-to-site hopping rate is controlled 
by quantum tunneling. In this temperature regime the 
de Broglie wavelength of a N adatom of mass m, i.e. 
Ab = 2nh/ \jZmk^T , is comparable to the distance sep- 
arating adjacent adsorption sites on the Au(lll) sur- 
face, therefore making quantum tunneling possible. For 
the fundamental energy level, the tunneling jump rate, 
Ttunneij can be derived for a one-dimensional parabolic 
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double-well potential and approximated by [33J 

rtunne^^^exp^-— j, (8) 

where uj = ^2E a /mb 2 and b is the barrier width [b = 
0.8 A for Au(lll)]. Therefore, the activation energy 
E a = 0.17 eV corresponds to r tunnc i = 66.9 Hz. If 
the calculated E a is considered to be accurate within 
±10 meV, r tU nnei can vary in the range ps 31.9—143.3 Hz. 

D. Work functions and surface dipoles 
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FIG. 3. Au(lll) slab models covered by a one-sided mono- 
layer (1 ML) of (a) He atoms and (b) N atoms. The varia- 
tion of the corresponding plane averaged electrostatic poten- 
tial V(z) along the z-axis normal to the surface is represented 
for (c) He-covered and (d) N-covered Au(lll) slabs, together 
with a clean Au(lll) slab reference. The Fermi energy Ey 
(horizontal dashed line) and the work functions (vertical ar- 
rows) of the clean and adsorbate-covered slabs, Wclean and 
Wads, respectively, are also represented. 

The work function, IT, is defined as the minimum en- 
ergy required to remove an electron from a solid to the 
vacuum region in the vicinity of the solid surface and is 
given by: 

W = V{oo) - E Fl (9) 

where V(oo) is the plane-averaged electrostatic potential 
in the vacuum at a distance where the microscopic po- 
tential has reached its asymptotic value and Ep is the 
Fermi energy. 

The electrostatic potential V(x, y, z) on a grid in real 
space can be obtained from a self-consistent electronic 
structure calculation using a plane wave basis set. As- 
suming that the surface normal is oriented along the z- 
axis, one can define a plane averaged potential 

V(z) = jJJ V(x,y,z)dxdy, (10) 



where A is the supercell surface area. The asymptotic 
value V(oo) can be extracted by plotting the variation 
of V as a function of z, as shown in Fig. [3] for a clean 
Au(lll) surface and for a Au(lll) surface covered by 
1 ML of He and N atoms. 

The calculated Fermi energy and electrostatic potential 
in the vacuum for the clean Au(lll) surface are Ep — 
—2.36 eV and V(oo) = +3.21 eV. This corresponds to a 
work function W c iean = 5.57 eV, in good agreement with 
the experimental value of 5.50 eV recently measured by 
Broker et al. for this Miller index plane [34] . 

Adsorption of one monolayer of He atoms at bridge 
sites and N atoms at fee sites onto the Au(lll) surface 
results in F(oo) = +3.24 eV and F(oo) = +4.68 eV, 
respectively. Thus, the work functions of He- and re- 
covered Au(lll) surfaces are W a( j s = 5.60 and 7.04 eV. 

In order to analyze the change of the work function 
upon atomic adsorption, we define the variation AIT = 
tTads — Wciean- For the weakly physisorbed He monolayer 
this variation is negligible (ATT = 0.03 eV), unlike in the 
case of the chemisorbed N monolayer (AIT = 1.47 eV). 
The variation of the work function results from the 
change in the surface electric dipole caused by adsorption 
of the adatoms. Simple electrostatics gives the relation 
[35] 

A„_3. <U> 

where A is the surface area taken up by one adatom, 
eo is the electric permittivity of free space and A/x is the 
change in surface dipole that occurs upon atomic adsorp- 
tion, normalized per adatom. A/i corresponds to the z- 
componcnt of the dipole moment directed along the sur- 
face normal, since only this component affects the work 
function. Since four adatoms form 1 ML covering the su- 
percell surface area, the surface area by adatom can be 
approximated by A sa A/ A. We can now estimate the 
induced dipole moment for the case of He and N adsor- 
bates. Using a DFT unit cell area of A ~ 41 A 2 we find 
^ Hc » 0.03 D while ^ N 1.60 D. 

It should be noted that the major contribution to the 
surface dipole results from the charge reordering associ- 
ated with the formation of the chemical bonds between 
the metal surface and the adatoms. This contribution is 
foremost determined by the nature of the chemical bonds, 
but can also be modified by the packing density of the 
adatoms. 



E. Phonon Density of States in Presence of the 
Monolayer Adsorbates 

Phonon density of states (PDOS) were calculated by 
solving the dynamical matrix for bulk Au, clean Au(lll) 
surface, and adsorbates (e.g., He or N) on the Au(lll) 
surface as shown in Fig. [4] A (2 x 2 x 2) supercell was 
adopted to obtain the force constant matrix of bulk Au 
that can be derived from the Hellmann-Feynman forces 
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obtained from the DFT calculations using VASP [IS] , 
The calculated PDOS of the bulk Au shows two main 
peaks represented by "T" and "L" that are in good agree- 
ment with previous experimental results |36j . According 
to Lynn et al. |36j . the longitudinal (L) and transverse 
(T) phonon modes of bulk Au are 4.61 THz and 2.75 
THz, respectively, at a high symmetry point (X) in the 
Brillouin zone. 
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FIG. 4. Phonon density of states (PDOS): (a) bulk Au, (b) 
bare Au(lll) surface, (c) 1ML He-covered Au(lll) surface, 
and (d) 1ML N-covered Au(lll) surface. The curves in black 
are the calculated total PDOS and the shaded areas in red, 
blue, and green represent the partial PDOS projected to the 
surface atoms, He atoms, and N atoms, respectively. 

The calculated surface PDOS using the (2 x 2) super- 
cells are depicted in Fig. [4jb) for a clean Au(lll) sur- 
face as well as for 1ML of adsorbates (He or N) on the 
Au(lll) surface. There are three prominent peaks in the 
calculated PDOS of the clean Au(lll) surface, mainly 
contributed by the partial PDOS projected onto the sur- 
face atoms (red shaded area) . The previous experimental 
study has identified four surface modes (2.31 THz, 3.5 
THz, 4.0 THz, and 4.3 THz) at a high symmetry point 
(K) in the surface Brillouin zone [37] . 

The calculated PDOS of the 1ML He-covered Au(lll) 
surface suggests a very weak interaction between He and 
the metal surface, providing no evidence of stretching or 
wagging modes of He atoms with respect to the metal 
surface. However, the partial PDOS projected onto the 
He atoms (blue shaded area) reveals the possible lattice 
modes of the 1ML He atoms physisorbed on the surface 
as shown in Fig. |4|c) . 

Contrasting with the 1ML He-covered Au(lll) surface, 
two additional peaks appear at high frequency above 5 
THz for the 1ML N-covered Au(lll) surface due to the 
wagging and stretching modes of N atoms attributed to 
the strong interaction with the metal surface, as shown in 
Fig.pjd). The green shaded area represents the partial 



PDOS projected onto the N atoms. 



IV. NOISE SPECTRUM WITH ML PRESENT 

Let us now study the impact of the two different types 
of MLs described above on the dipole fluctuation spec- 
trum of adatoms as discussed in Sec. [Til We first note 
that due to the high reactivity and low mass of H ad- 
sorbates - chosen in the previous section to reduce the 
DFT computational cost - the depth of the binding po- 
tential Uo,tare — 0.068 Hartree and vibrational frequen- 
cies ^io, bare/27r « 40 THz are high. Therefore, for H 
adatoms thermally activated processes at room temper- 
ature do not play a significant role. Instead we use the 
potentials shown in Fig. [2] together with the mass scal- 
ing relations in Eq. ([6| to evaluate the noise spectrum 
for a more realistic set of adsorbates with masses around 
m ~ 100 amu. 
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FIG. 5. The dipole fluctuation spectrum Syi(w) as a function 
of lo/To : au for He monolayer (black), N monolayer (red) and 
bare system (blue). For each system the solid curve corre- 
sponds to T = 50 K while the dashed curve corresponds to 
T = 150 K. The frequency ui is scaled by Fo.Au, the char- 
acteristic rate corresponding to He monolayer coverage. The 
spectrum is given in units of /iq/To,Au, where fio is the in- 
duced dipole moment for the ground state for each system. 

In Fig. [5] we plot the dipole fluctuation spectrum S^^uj) 
for a bare Au surface and in the presence of a He and 
N monolayer. For the results shown in Fig. [5] we use 
the potentials shown in Fig. [2] with an adatom mass of 
100 amu and two different temperatures of T = 50 K 
(solid) and T = 150 K (dashed). For the bare Au surface 
where z w 1.59 A these values correspond to v w /2-k rs 
4.5 THz, T^ Au /2tt m 2.2 THz and ratios k B T/hu w 
0.20 and 0.70 respectively. For a qualitative discussion of 
the strong modifications of the spectrum in the presence 
of a ML, we consider low temperatures, where the adatom 
potential can be approximated by only two vibrational 
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level system and 



SUMMARY AND CONCLUSIONS 



s n(u) = (Mi ^ Mo)" 



2r 



r 2 

1 



-hv 10 /k B T 



(12) 



In Sec. we have found that due to its low reactiv- 
ity, the He monolayer results in a significantly shallower 
well depth, Uo.He = 0.00055 Hartree, shifting the min- 
imum to zq^hc ~ 5 A. Both effects lower the vibra- 
tional frequency, vio t He ~ 0.4 THz, and lead to a dras- 
tic reduction of the characteristic phonon transition rate 
r ,_f/ e /27r sa 140 MHz. From Eq. ([12]) this results in an 



increase of the low frequency noise level, but reduces the 
noise at frequencies uj ^> Fq. 

The more reactive N monolayer results in a deeper po- 
tential well with Uo.n — 0.141 Hartree, while only slightly 
affecting the equilibrium distance, zo.w ~ 2 A. This 
leads to exactly the opposite effect, increasing the vibra- 
tional frequency, viq^n ~ 5.3 THz as well as Tq^n/2tt « 
3.9 THz. As also seen in Fig. [5j the stronger binding 
therefore reduces the noise in the low frequency regime. 
Additionally the low frequency regime of the spectrum 
extends over a much larger frequency range. 

Finally, let us briefly comment on the modified PDOS. 
Our estimates of r so far have been based on Eq. ffify, 
which assumes g(u)) oc oj 2 . This scaling ignores the de- 
tailed structure of the PDOS at high frequencies and 
overestimates Tq for V\q above ~ 1 THz. Although the 
variations in PDOS shown in Fig. [4] are in general less 
relevant and only important for adatoms with high vi- 
brational frequencies, it seems feasible to identify ML 
species, for which the PDOS is peaked at the vibrational 
frequencies v\q, thereby increasing Tq and suppressing 
the low frequency noise level. 



In summary, we have calculated the noise induced 
heating in ion traps due to randomly-fluctuating adatom 
dipoles in the presence of a single ML of atomic species on 
Au(lll) surface. Precise DFT calculations of surface po- 
tentials for physisorbed and chemisorbed ML species pro- 
vided us with accurate data for surface potentials, from 
which the effects on the noise could be evaluated. Wc 
have found that within our noise model, the presence of 
surface contamination can lead to opposite effects of en- 
hancing or reducing the noise level, depending on the re- 
activity of the ML species as well as the frequency range. 

Although an exact quantitative comparison between 
experiments, for example, in Ref. |12j and |13j . is beyond 
the scope of this work, our current analysis points the di- 
rection for a more refined understanding of anomalous 
heating in ion traps, requiring the combined knowledge 
of atomic surface physics as well as the modeling of differ- 
ent noise processes. In particular, the combined data on 
surface potentials, PDOS and diffusion rates obtained in 
this work could in the future serve as a common input to 
evaluate and compare different alternative noise models. 
This can lead to better understanding of the noise mech- 
anism, allowing ion traps to be used as exquisite probes 
of surface reactivity. 
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